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Identification of beyond-standard-models including WIMP dark matter is studied in 
four particle final state with a W boson pair and a WIMP pair at the International 
Linear Collider. Models with different spin structures give distinguishable production 
angle distributions. After the mass determination in each model, the production angle 
is reconstructed using the four momentum of W bosons with a back-to-back constraint. 
Three models of Inert Higgs, Supersymmetry and Little Higgs are considered. Discrim- 
ination power at 200 fb and 40 fb signal cross section with 500 fb -1 luminosity at 
y/s = 500 GeV is obtained. 
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H ' 1 Introduction 



A weakly interacting massive particle (WIMP) is a leading candidate for dark matter. Many 
models including particles of various spin arc proposed in order to explain existence of the 
WIMP dark matter. Using angular distribution, those models may be separated in collider 
experiments. 

In this study, we investigate the separation possibility in the process of e + e~ — > x + X~ ~^ 
X°X°W + W _ at the International Linear Collider (ILC) with full Monte Carlo simulation. \° 
and x ± denote the WIMP and the lightest charged particle (LCP) in the model, respectively. 
Center-of-mass energy of ILC is set to be 500 GeV. Particles in the model other than the 
LCP and the WIMP are assumed to be inaccessible at this energy. 

We consider three WIMP model candidates as shown in Table [TJ Since the three models 
have different spins of their particles (Inert Higgs (IH)[5]: scalar, Supersymmetry (SUSY) [3]: 
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Model 


Particle 


Name 


Spin 


Mass 


a at 500 GcV 


Inert Higgs 


LCP 
WIMP 









231.57 GcV 
44.03 GeV 


3.51 fb 


Supersymmetry 


LCP 

WIMP 


x ± 

X? 


1/2 
1/2 


231.57 GeV 
44.03 GeV 


384 fb 


Little Higgs 


LCP 

WIMP 




1 
1 


231.57 GeV 
44.03 GeV 


364 fb 



Table 1: Models and cross sections used for this study. In the analysis Cross sections are 
normalized to 40/200 fb so the cross section listed in this table is not used. 



fermion, Little Higgs (LH)[S]: vector), they are good samples for the model identification 
study. Mass is set to be the same as m(x ) = 231.57 GeV and m(x°) = 44.03 GeV for each 
model by tuning its model parameters. Cross section is normalized to 40 fb and 200 fb for 
the analysis. 



2 Simulation and Reconstruction 

The event generation is performed using Physsim[2] event generator for events of three 
sample models and Whizard[5] generator for the Standard Model (SM) background sample. 
The initial-state radiation and beamstrahlung are included in the event generations. The 
RDR[7] nominal beam parameters are used while the finite crossing angle is ignored. We 
assume no beam polarization. 

According to the assumed cross section of 200 fb, 0.1 million signal events are generated 
for each model. The SM background sample is prepared for International Large Detector 
(ILD) Letter of Intent (Lol) jS]. Its statistics corresponds to 0.1 - 500 fb -1 luminosity 
(depends on processes) with about 12 million events in total. 

ILD full simulation code (Mokka[H] and MarlinReco[Tn]) is used for the full Monte- 
Carlo (MC) simulation and event reconstruction. ILD_00 detector geometry, which is the 
standard geometry for ILD Lol is used for the simulation. The geometry includes a time 
projection chamber with silicon devices for tracking and vertexing, and highly granular 
electromagnetic and hadronic calorimeters for particle flow calorimetry along with 3.5 Tesla 
magnetic field. The simulation code includes many realistic features, gaps in the sensitive 
regions and estimates of dead material due to cables, mechanical support, cooling and so 
on. The central part of the reconstruction is the particle flow algorithm PandoraPFA [IT] , 
which forms charged and neutral particle candidates (particle flow objects: PFOs) from 
tracks and calorimeter clusters. The resulting list of PFOs for each event is forced into a 4- 
jet configuration using the Durham algorithm |12|. Neural- net based flavor tagging algorithm 
LCFIVertex[T3] is applied for the jets after the clustering and vertex finding with ZVTOPpB] 
algorithm. As a final step of the reconstruction, a constrained kinematic fit |15j . which 
requires the two dijet masses of the event to be equal, is performed on each event. All three 
possible jet pairings are tested and the pairing with least x 2 value for the kinematic fit is 
adopted. 
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3 Signal Selection 

We use signal events with both W decaying into two quarks (qqqq events), which has about 
46% branching ratio, since W energy must be fully reconstructed for the mass determination 
and the production angle reconstruction. The signature of the fully hadronic decay from our 
target processes is '4-quarks + missing'. The main SM background processes are W- and 
Z-pair production with fully hadronic decay, top-pair production with one leptonic decay of 
W, 77 — > WW processes, WW ' Z process with the Z decays to a impair, etc. Scmi-lcptonic 
decay of the target processes is also background for our study. 

The following selection is applied to all samples to reject the major part of SM and 
semi-leptonic decay background: (i) Number of tracks should be larger than 20 and each jet 
has to contain at least two tracks in order to eliminate pure leptonic events, (ii) The visible 
energy of the event E v - ls should be between 80 and 400 GeV which can remove most of 2- 
photon and 2/4/6 quark events, (iii) each jet should have a reconstructed energy of at least 
5 GeV and a polar angle 6 fulfilling |cos(0j ct )| < 0.99 to ensure proper jet reconstruction, 
(iv) Requiring the distance parameter of the Durham jet algorithm|12| for which the event 
changes from 4-jet to 3-jet configuration, j/34 to be larger than 0.001 in order to reject most 
of 2-jet events, (v) No lepton candidate with an energy larger than 25 GeV is allowed in 
order to suppress semi-leptonic events, (vi) |cos(#)| of the missing momentum should be 
smaller than 0.9 and |cos(0)| summed up for all jets should be smaller than 2.6 in order 
to eliminate most of SM events which are concentrated in the forward region, (vii) B-tag 
probability from LCFIVertex output summed up for all jets should be smaller than 1 to 
remove events with b quarks, (viii) The kinematic fit constraining the two dijet masses to 
be equal should converge for at least one jet pairing to ensure integrity of the fit result, (ix) 
The di-jet mass obtained by the kinematic fit should be between 65 and 95 GeV to select 
two-W events. 
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Figure 1: Cut plots in 200 fb signal cross section. The labels (ii), (vii), (ix) correspond to 
the cuts described in the text with the same labels. 

The obtained reduction due to these cuts is shown in Table [2] and some of the cut plots 
are shown in Figure Q] Clear peaks at di-jet masses of 80 GeV can be seen in the signal 
distributions of Figure [1] (ix), which are from two W bosons. Acceptances of signal events 
are 59.8, 58.3 and 58.9 % for IH, SUSY and LH models, respectively. Signal purities after 
the cuts are 78.1, 77.3 and 77.6 % in 200 fb cross section and 42.1, 40.9 and 41.6 % in 40 fb 
cross section, respectively. 
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Processes 


No cut 


before mass cut 


all cuts 


Signal 


Inert Higgs (hadronic decay) 


46746 


33422 


27943 


SUSY (hadronic decay) 


45539 


31942 


26552 


Little Higgs (hadronic decay) 


46466 


32681 


27361 


Model bkg. 


Inert Higgs (other decay) 


53252 


1973 


185 


SUSY (other decay) 


54462 


1930 


139 


Little Higgs (other decay) 


53532 


1770 


239 


SM bkg. 


qqqq (WW, ZZ) 


1.88e+06 


8360 


3223 


qqHv [WW) 


2.35c+06 


6796 


1885 


qqqqlv (tt) 


125204 


2754 


626 


qqivlv 


31630 


103 


39 


11 -> qqqq 


26356 


666 


523 


qqqqvv (WW Z) 


4158 


1843 


681 


qqvv (ZZ) 


117808 


8547 


134 


qq 


6.29c+06 


1518 


373 


11 -» qq 


7.97c+06 


914 


58 


Other background 


3.42c+09 


727 


108 



Table 2: Event number before and after selection cuts for 200 fb production cross section, 
normalized to 500 fb -1 and no beam polarization. 



4 Mass Determination 



The masses of new particles can be obtained via the energy spectrum of the W boson 
candidates. The energy of W bosons by production of new particles have upper and lower 
kinematic limits from which the masses of the new particles can be derived. Figure [2^a) 
shows the W energy spectrum for each model with SM background. Mass edges can be seen 
in the distribution of every model. 

The edge positions of the energy spectra are obtained by fitting. The fitting function is: 
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Figure 2: (a)W energy distribution of signal and background. (b),(c) Results of the mass 
fit for 200 fb and 40 fb signal cross section after background subtraction, respectively. 
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Model 


a 


Fitted left edge 


Fitted right edge 


LCP 8m 


WIMP 6m 


True edges 




96.32 GeV 


174.78 GeV 






Inert Higgs 


200 fb 

40 fb 


97.00 ± 0.13 GeV 
97.86 ± 0.38 GeV 


174.22 ± 0.19 GeV 
176.22 ± 0.84 GeV 


0.12 GeV 
0.39 GeV 


0.57 GeV 
2.23 GeV 


Supersymmetry 


200 fb 
40 fb 


97.12 ± 0.15 GeV 
96.10 ± 0.54 GeV 


173.46 ± 0.21 GeV 
176.02 ± 0.85 GeV 


0.14 GeV 
0.50 GeV 


0.64 GeV 
2.50 GeV 


Little Higgs 


200 fb 

40 fb 


96.94 ± 0.16 GeV 
99.56 ± 0.48 GeV 


174.39 ± 0.10 GeV 
174.67 ± 0.34 GeV 


0.13 GeV 
0.41 GeV 


0.50 GeV 
1.56 GeV 



Table 3: Result of the mass fit. Refer to the text for 8m calculation. 



/(x;io-i,&o-2,0o-i,r) = / (b 2 t 2 + b 1 t + b Q )V(t- x, atf + <r a ,T)dt 



(1) 



where V(x, a, T) denotes the Voigt function of Gaussian width of a and Lorenzian width 
of r. After the fit with all parameters free (8 free parameters), we fixed all the parameters 
except to an d t\ (which gives the edge positions) and refit to obtain shape-independent mass 
resolution. Figure [2jb)(c) gives the fitting result with 200 fb and 40 fb signal cross section, 
respectively. The fitting result is summarized in Table [3] While the center value of the 
fitting result deviates from the expected true edge position, it can be corrected using Monte 
Carlo samples in the real experiment. 

The masses of LCP and WIMP can be calculated from the edge positions. Error of the 
edge positions can be translated to the error of the masses 8m. For the calculation, true edge 
positions are used instead of the (deviated) center value of the fitting result. The obtained 
mass error is also shown in Table G3 

Since we can see clear difference among three models in the W energy distribution, the 
model separation may be possible with the difference. However, the difference is considered 
to come from the specific model and not from the general spin structure, we do not use the 
difference to identify the model in this study. 

5 Angular Distribution for LCP Pair Production 

Separation of the models is possible by comparing distributions of production angles of the 
LCPs. To derive the production angles, a quadratic equation is solved with the masses of new 
particles and the momenta of the W bosons with the assumption of back-to-back ejection 
of two LCPs. The equation gives either two solutions which contain one correct production 
angle or no solutions when the discriminant of the equation is negative. The unphysical 
negative discriminant comes from misreconstructing W momenta or imperfect back-to-back 
condition of the two LCPs mainly due to initial state radiation. 23.9% (IH), 20.8% (SUSY), 
23.7% (LH) and 64.4% (SM background) of the events have negative discriminant and are 
cut off for the following analysis. 

Figure [3] shows production angle distributions. 1-dimensional results show the visible 
difference between three models that IH events concentrate in the central region while SUSY 
events are almost flatly distributed. LH events give intermediate profile. 2-dimensional 
results are actually used to calculate the separation power. To quantify the difference, we 
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Production angle, MC truth 
2500 




(a) 



(b) 







- 120 — 




a 



(d) 



(e) 



(f) 



Figure 3: Production angle distributions with 200 fb cross section, (a) and (b) show MC 
and reconstructed 1-dimcnsional distribution. Both of two solutions are included, (c)-(f) 
give 2-dimcnsional distribution of IH, SUSY, LH and SM background, respectively. 



Data / Template 


IH 


SUSY 


LH 


Inert Higgs 


- 


86 


20 


Supersymmetry 


71 


- 


31 


Little Higgs 


23 


27 


- 



D / T 


IH 


SUSY 


LH 


IH 


- 


11.7 


2.6 


SUSY 


9.9 


- 


4.5 


LH 


2.9 


4.1 


- 



(a) 200 fb 



(b) 40 fb 



Table 4: Separation power P among three models with the 2-dimensional production angle 
distribution. The numbers include uncertainty of ~ 2 for 200 fb and ~ 0.3 for 40 fb signal 
cross section due to limited MC statistics. 



compare the production angle distribution for one model (dubbed as "data" ) against another 
model ( "template" ) . The chi-square variable x 2 an d the separation power P is given by: 



bins (n - T) 2 - (t 2 



MC,i 



Di 



P 



N 



(2) 



where Di and TJ are number of data and template events (including SM background) in the 
ith bin, <7mc,i is the standard deviation of MC statistics in both data and template samples, 
and N is the number of bins (210). The template events are normalized to the integral of 
the data events before calculating \ 2 . 

Table |4] shows the obtained separation power P. (a) shows enough separation power 
> 20<7 should be obtained with 200 fb signal cross section. In 40 fb case (b), 2.6 - 2.9 a 
separation power is expected between IH and LH models while > 4cr separation is performed 
in other pairs. 

Since the production angle reconstruction uses the masses of the new particles, deviation 
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Fit with the power of 0.5 
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\fs [GeV] 
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(a) 



(b) 



Fit with the power of 1 .5 



£ 10 




10 20 30 40 50 
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Figure 4: (a) Dependence of the cross section on the beam energy, normalized to 40 fb at 
y/s = 500 GeV. Error bars are given assuming 50 fb -1 data at each point, (b)-(c) Result of 
the x 2 nt. (b) (y/s- y/^) 1/2 case, (c) (y/s ■ 



so) 3 / 2 case. 



of mass determination can cause degradation to the production angle study. To investigate 
systematic effect from the mass uncertainty, we produce more template samples with shifted 
MC masses of LCP and WIMP. The amount of mass shift is determined as the same as the 
mass measurement error at 40 fb signal cross section (shown in Table [3]). Both negative and 
positive shifts are considered (four more templates per model are produced) . The separation 
power obtained by the shifted templates is consistent with the result of the original templates 
within the MC statistics. From that, we can conclude that the effect of the mass uncertainty 
to the model separation is no more than the uncertainty of the MC statistics (~ 2 in 200 fb 
and - 0.3 in 40 fb). 

6 Threshold Scan 



Another strategy to identify the models is to change center-of-mass energy and monitor 

difference of the cross section. Figure 0] (a) shows the cross section dependence of each 

l li 
model on Eqm ■ A clear difference can be seen as the SUSY model of E^ M dependence and 

the other two of E C ' M dependence. 

To obtain separation power with the threshold scan, we perform a toy-MC study in which 
the measured cross section is fluctuated using signal and background statistics obtained 
from the full-MC study. The cut efficiency and background cross section are assumed to be 
identical to the 500 GeV case in all y/s . We assume a 3-point scan of yfs =470, 500 and 
530 GeV with 50 fb -1 integrated luminosity each. Cross section is set to 40 fb at 500 GeV. 

For the separation, we calculate the \ 2 value of the fit of, 



cr(\/s, n) = a( y/s - y/loj 7 



= 1/2,3/2 



(3) 



where a and y/so are the free parameters for each model. Figure 2] (b)-(c) shows the x 2 
distributions. With the power of 1/2 fit (b), clear separation is obtained between SUSY and 
other two models. For example, 99.4% of the SUSY events are within x 2 < 12 while 0.11 
and 0.38% of the IH and LH events remain for the same \ 2 region. The power of 3/2 fit (c) 
has less separation power, and separation between IH and LH is almost impossible by the 
threshold scan. 
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7 Summary 

We study possibility of the identification of WIMP models which have different spin struc- 
ture with W + dark matter pair production at s/s = 500 GcV ILC. With 200 fb cross 
section, separation power among three models is strong enough with the production angle 
distribution and ~ 3er separation can still be obtained with 40 fb cross section. We also show 
that the threshold scan is a very powerful method for lower statistics, while in some cases 
(ex. IH and LH) behavior around the threshold is similar and separation is not practical. 
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